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The band gaps, longitudinal and transverse effective masses, and deformation potentials of ScN 
in the rock-salt structure have been calculated employing GoWo-quasiparticle calculations using 
exact-exchange Kohn-Sham density functional theory one-particle wavefunctions and energies as 
input. Our quasiparticle gaps support recent experimental observations that ScN has a much lower 
indirect band gap than previously thought. The results are analyzed in terms of the influence of 
different approximations for exchange and correlation taken in the computational approach on the 
electronic structure of ScN. 
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'. I. INTRODUCTION 

Scandium nitride (ScN) is emerging as a versatile ma- 
terial for promising technological applications. As part 
of the transition metal nitride family it initially gener- 

[ ated interest for potential applications as wear resistant 
and optical coatings due to its mechanical strength, high 

■ melting point of 2600° C^, and high hardness {H = 21 
GPa) with respect to load deformations^. ScN crystal- 
lizes in the rock-salt phase with a lattice parameter of 
4.50 A^. The octahedral bonding arrangement provides 
a much more favorable environment for the incorporation 
of transition metal atoms like Mn or Cr than the tetrahe- 
drally coordinated III-V semiconductors, which have up 

' until now been popular candidates for spintronic materi- 
als. Successful incorporation of Mn into ScN has been 
demonstrated^ and ab initio calculations predict Mn- 
doped ScN to be a dilute ferromagnetic semiconductor^. 
Moreover, ScN has a lattice mismatch of less than 2% to 

i cubic gallium nitride (GaN). This makes ScN structurally 
compatible with the group-IIIA nitrides^'i'^'^iiS'iiii^!^ - 
an important technological material class, in particular 
for applications in optoelectronic devices. Alloying ScN 
with GaN^ii^iiiii^ might provide a viable alternative to 
InGaN alloys for use in light emitting devices or solar 
cells. In addition multifunctional devices are conceiv- 
able if the strong electromechanical response predicted 
for hexagonal ScNi^ can be utilized. 

The electronic band structure of ScN - a key quan- 
tity for the design of optoelectronic devices - has been 
difficult to access both experimentally and theoretically. 
Early experiments were hampered by various complica- 
tions in growing films with well defined crystalline ori- 
entation, stoichiometry, low background carrier concen- 
tration, and surface roughness. For a detailed discussion 



we refer to e.g. Ref. [H. Recent advances in growth 
techniques have led to a systematic improvement of the 
material's quality^*'. Employing optical spectroscopy and 
photoemission. Gall et alr^ concluded that ScN is a semi- 
conductor with an indirect T — X band gap {Eg~'^) of 
1.3±0.3eV. The sizable error bar of 0.3 eV has been 
mainly attributed to the large background carrier concen- 
tration of ~ 5 X 10^" cm^"^ causing an apparent increase of 
the band gap due to the Burnstein-Moss shift^^. Reduc- 
ing the electron carrier concentration to 4.8xl0^^cm~^ 
and combining tunneling spectroscopy and optical ab- 
sorption measurements, Al-Brithen et ali^ were able to 
reduce the error bar and found a value for Eg~-^ of 
0.9±0.1eV. 

Early Kohn-Sham density functional theory (KS- 
DFT) calculations employing the local-density (LDA) or 
Xa approximations predicted ScN to be a semimetal 
with a small negative band gap between -0.01 and 
-0.21 eV— In order to overcome the well known 
underestimation of the LDA band gap, more ad- 
vanced exact-exchange (OEPx(cLDA))i^ and screened 
exchangers, calculations have been performed, and 
showed that ScN is a semiconductor with an indirect F to 
X band gap, in accord with experimental evidencei^i^. 
However, the calculated band gap of 1.60 eV found in 
both studies is significantly larger than the most recent 
experimental value of 0.9±0.1eV— . 

In order to shed light on this discrepancy we have per- 
formed quasiparticle energy calculations in Hedin's GW 
approximation^^, which is a well established technique to 
calculate accurate band structure energies and currently 
the choice for computing quasiparticle band structures of 
solid o^^i^^'^^ . The quasiparticle calculations predict ScN 
in the rock-salt phase to have an indirect band gap be- 
tween the F and X point of 0.99±0.15eV, strongly sup- 
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porting recent experimental findings. In addition we have 
also determined the direct band gaps and other electronic 
structure parameters relevant for device simulations: the 
volume deformation potentials of the main band gaps 
and the longitudinal and transverse effective masses of 
the conduction band at the X point. The effective mass 
has previously been calculated at the level of the LDASl, 
but to the best of our knowledge only one experimental 
study has reported a conduction band effective mass for 
ScN (between 0.1 and 0.2 TOq)^ so far. For the defor- 
mation potentials no experimental or theoretical data is 
available, yet. 

Most commonly, the Greens function Gq and the 
screened potential Wq required in the GW approach 
(henceforth denoted GqWq) are calculated from a set 
of KS-DFT single particle energies and wave functions 
{ei, (pi}. Since Go and Wo are not usually updated by the 
quasiparticle wave functions and energies {ef^, 'pf^} 
a self-consistent manner, the quasiparticle energies de- 
pend on the approximation used to calculate the input 
datc^i^s^i^. 

Originally, GqWq calculations were based on LDA 
data (LDA-GqWo) and were found to accurately pre- 
dict band gaps of sp-bonded semiconductors (with a typ- 
ical error bar of ^O-leV)^"*. However, complications 
arise when the LDA-GqWo approach is used to calcu- 
late the electronic structure of semiconductors with neg- 
ative LDA band gaps'^^'^'^ or when occupied shallow semi- 
core d bands are treated as valence in the pseudopoten- 
tial framework^^i^i^^. For such semiconductors, GqWq 
calculations based on OEPx(cLDA) data (OEPx(cLDA)- 
GqWq) have been found to provide a reliable tool to ob- 
tain band gaps with an accuracy of 0.1 eV^^'^^. 

The key to the improved description in the 
OEPx(cLDA)-GoWo approach can be found in the treat- 
ment of exchange. In the exact-exchange KS approach 
the formal expression for the total energy is the same as 
in Hartree-Fock. The difference between the two meth- 
ods lies in the potential felt by the electrons: in Hartree- 
Fock the exchange potential is non-local, whereas in the 
exact-exchange KS approach it is local and constructed 
to be the variationally best local ground state poten- 
tial to the non-local Hartree-Fock exchange potentiai^. 
Like Hartree-Fock, the exact-exchange KS approach is 
therefore free of self-interaction, but since the eigenval- 
ues are solutions to a local potential they are in general 
closer to (inverse) photoemission data for semiconduc- 
tors than the Hartree-Fock single particle energie o^^i'^^'^" . 
The GW formalism, on the other hand, goes beyond 
the KS approach and describes the interaction of weakly 
correlated quasiparticles by means of a non-local, en- 
ergy dependent self-energy. It takes the form of the 
non-local exchange potential encountered in the Hartree- 
Fock approach, which is screened by correlation in the 
random-phase approximation (RPA). To elucidate the ef- 
fects of exchange and correlation on the electronic struc- 
ture of ScN we therefore first analyze the influence of 
exchange by comparing LDA (GGA) and OEPx(cLDA) 



calculations, before turning to the difference between 
OEPx(cLDA) and GoWq. 

The paper is organized as follows. Sec. [H] describes 
our computational approach. The results are presented 
and discussed in Sec. IIVI Finally, a summary is given in 
SecHVl 



II. COMPUTATIONAL METHOD 

The KS-DFT calculations have been performed with 
the ab initio pseudopotential (PP) plane- wave code 
SPHIngX**^. A consistent set of norm-conserving scalar- 
relativistic PPs has been used for each of the exchange- 
correlation functional (LDA, GGA and OEPx(cLDA)). 
The OEPx(cLDA)-PPs have been constructed as de- 
scribed in Ref. 42 and the LDA and GGA ones using 
the FHI98PP code"^. All PPs for both Sc and N have 
been generated according to the TrouUier-Martins opti- 
mization scheme'*'* and have then been transformed into 
the separable Kleinman-Bylander form^^. Unlike in the 
group-IIIA-nitrides GaN and InN the cation 3s and 3p 
states have moderate binding energies in ScN and the Sc 
3p and to a lesser degree also the Sc 3s electrons couple to 
the upper valence bands as a partial charge density anal- 
ysis reveals (see Section IlII Al and also Fig. 4 in Ref. [23 ). 
For Sc, the entire semicore shell (3s, 3p and 3c? states) 
has therefore been treated as valence with an ionic con- 
figuration [Ne]3s23/3d*. Following Ref. \M, only the 2s 
and 2p components have been included for N. Adding a 
d component for N yields negligible effects on the calcu- 
lated band structure and total energies of ScN. For Sc, 
we have chosen a core radius of 1.4, 1.4, and 1.8 bohr 
for the s, p and d orbitals, respectively. For N, a com- 
mon core radius of 1.5 bohr has been adopted. The s (p) 
component is taken as the local component for Sc (N). 
Only one projector per angular momentum channel has 
been used, i.e. 3p and 3d for Sc. We have verified that 
this procedure does not compromise higher lying states 
in the same channel by ensuring that the eigenvalues of 
the 4s and 4p levels in the pseudo-atom reproduce those 
of the all-electron calculation. These pseudopotentials 
have been carefully tested (see also below) and are free 
of ghost states^. 

For the LDA calculations we have used the Ceperley- 
Aldcr'** exchange-correlation data as parametrized by 
Perdew and Zunger"*^. The GGA calculations have been 
performed with the Perdew-Burke-Ernzerhof^° func- 
tional. In the OEPx(cLDA) calculations, the exchange 
energy and potential have been treated exactly and cor- 
relation has been added on the LDA level. Throughout 
the paper, the combination of exact-exchange and LDA 
correlation will be referred to as OEPx(cLDA). 

The Go Wo calculations have been performed em- 
ploying the GW space-time approach^, in the gwst 
implementation^i^i^. The Kohn-Sham eigenvalues and 
wavefunctions (e^ and in either OEPx(cLDA) or LDA 
are used as input to construct Go and Wq. Head and 
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Approach 


ao (A) 


Bo (GPa) 


B'o 


Present work 


PP-PW(LDA) 


4.455 


221 


4.27 


PP-PW(GGA) 


4.533 


196 


4.36 


Other theoretical calculations 


FP-LAPW(LDA)^^ 


4.42 


235 




FP-LAPW(LDA>2i 


4.44 


220 




FP-LAPW(GGA)S 


4.50 


201 




FP-LAPW(GGA)'^i 


4.54 


201 


3.31 


Experiment^ 


4.501 


182±40 





TABLE I: Structural parameters of ScN: lattice constant (ao), 
bulk modulus (Bo) and its derivative {B'q) calculated in our 
pseudopotential, plane-wave (PP-PW) approach compared to 
previous results obtained with the all-electron full-potential 
linearised augmented plane-wave approach (FP-LAPW) and 
experiment. 



wings of the dielectric matrix (which converge slowly with 
respect to the k-mesh) have been calculated using a fine 
10 X 10 X 10 Monkhorst-Pack (MP) mesb^^. We find that 
using an offset of [5,5,5] yields faster convergence with 
respect to the number of fc-points, because the fc-point 
set then contains fewer high-symmetry points. Contri- 
butions arising from the non-local part of the pseudopo- 
tential are fully taken into account^. A regular 4x4x4 
MP mesh centered on the F-point then proves to be suf- 
ficient for the full GqWq calculations. 

Brillouin zone integrations in the DFT calculations 
have been performed on a 4 x 4 x 4 MP mesh. In all cal- 
culations an energy cutoff of 80 Ryd is used for the plane 
wave expansion of the wavefunctions. For the indepen- 
dent particle polarizability xo^ which enters in the calcu- 
lation of the OEPx(cLDA) potential (see e.g. Ref. [soh . 
an energy cutoff of 55 Ryd gives converged results. Con- 
duction band states up to the same energy cutoff have 
been included in the calculation of the electronic Green's 
function in the OEPx(cLDA) as well as in the GW calcu- 
lations. These parameters yield converged KS and quasi- 
particle energies to within 0.05 eV. 

We have tested the reliability of our pseudopotential 
calculations against all-electron calculations in the full 
potential linearized augmented plane wave (FP-LAPW) 
approach, where available. For ground state calcula- 
tions on the level of LDA and GGA we find good agree- 
ment for the structural properties, as Table U illustrates. 
LDA slightly underestimates the lattice constant of ScN 
compared to experiment^ by approximately 1%, whereas 
GGA results in a slight overestimation. 

The issue of performing self-consistent GW calcula- 
tions is still a matter of debate^iSi^i^i^. Unlike in 
DFT, a self-consistent solution of the full set of equations 
for the self-energy in many-body perturbation theory 
would go beyond the GW approximation and successively 
introduce higher order electron-electron interactions with 
every iteration step. Solving the GW equations self- 
consistently is therefore inconsistent if no higher order 



Approach 




El-- 


Eg 


Present work 


OEPx(cLDA)-GoVKo 


3.51 


0.84 


1.98 


LDA-GoW^o 


3.71 


1.14 


2.06 


[Go Wo] average 


3.62 


0.99 


2.02 


OEPx(cLDA) 


4.53 


1.70 


2.59 


GGA 


2.43 


-0.03 


0.87 


LDA 


2.34 


-0.15 


0.75 


Other theoretical work 


OEPx(cLDA)^'^ 


4.70 


1.60 


2.90 


sX^ 




1.58 


2.41 


Experiment 


Ref. 


~3.8 


1.30 


2.40 


Ref. 18 




0.9±0.1 


2.15 









TABLE IL Calculated and experimental band gaps (Eg) of 
ScN (in eV). sX denotes previous screened exchange calcula- 
tions and [Go Wo] average the arithmetic average between the 
OEPx(cLDA)-GoWo and LDA-GoWq results (see text). 



electron-electron interactions are included. It was first 
observed for the homogeneous electron gas^^ that the 
spectral features broaden with increasing number of iter- 
ations in the self-consistency cycle. Similarly, for closed 
shell atoms the good agreement with experiment for the 
ionization energy after the first iteration is lost upon 
iterating the equations to self-consistency^. Imposing 
self-consistency in an approximate fashion^i^i^!^ is 
not unique and different methods yield different results. 
Since the issue of self-consistency within GW is still dis- 
cussed controversially, we refrain from any self-consistent 
treatment and remain with the zeroth order in the self- 
energy (GqWo). We argue (see Section ITlI Al) that in the 
case of ScN the error bar resulting from this approxima- 
tion is only of the order of 0.15 eV. 



III. RESULTS AND DISCUSSION 
A. Electronic band structure 

The quasiparticle band structure of ScN is calcu- 
lated employing both the OEPx(cLDA)-GoWo and LDA- 
GqWq approach. To understand the effect of the start- 
ing point on the GqWq calculations we first analyze 
the KS band structure using three levels of approxi- 
mations for the XC potential (OEPx(cLDA), GGA and 
LDA). For a meaningful comparison between the results 
of these calculations among themselves and with exper- 
iment, these electronic structure calculations have been 
performed at the experimental equilibrium volume. As 
an example, we show in Fig. [1] the electronic band struc- 
tures of ScN in LDA, OEPx(cLDA) and OEPx(cLDA)- 
GqWq. Table HIl summarizes the calculated band gaps 
{E^^^, Ef-^ and E^-^), previous 0EPx(cLDA)i5 and 
screened exchange^ results and experimental dataF^i^^. 
For the following discussion we consider only the latest 
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FIG. 2: Schematic diagram of the pd bonding in rock-salt ScN: 
the Sc Sd states are split into t2g and eg states by the crystal 
field. The t2g states then form bonding and anti-bonding 
bands with the N 2p states and the eg states non-bonding 
bands that lie higher in energy. The occupation of the relevant 
states and bands is shown by circles. 
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FIG. 1: Electronic band structure of rock-salt ScN in LDA, 
OEPx(cLDA) and OEPx(cLDA)-GoW^o. The top panel shows 
the upper valence and lower conduction bands, aligned at the 
top of the valence band. The bands in the lower panel have 
mainly N 2s (around -13 eV), Sc 3p (around -30 eV) and Sc 
3s character (around -50 eV). 



experimental data of Al-Brithen et al. on low background 
carrier samples^^ as a reference. 

In ScN the scandium atom donates its two 4s and sin- 
gle 3rf electron to the nitrogen atom. According to the 
bonding analysis of Harrison and Straub^, the five d 
states of Sc hybridize with the three valence p states 
of the neighboring N atoms in the rock-salt structure 
of ScN, forming three p-like bonding, three d-like anti- 
bonding t2g and two rf-like non-bonding Cg bands. The 
bonding scheme together with the electron filling of these 
bands is sketched in Fig. [2l Performing a partial charge 
density analysis we have confirmed that the upper three 
valence bands in the DFT calculations correspond to the 
bonding states and originate mainly from the N 2p states 
with some admixture of the Sc 3d states, while the low- 
est conduction band are the anti-bonding t2g states with 
Sc 3c? character. The two bands derived from the non- 
bonding Cg states are around 1.2 eV higher in energy. 
This assignment is consistent with the partial density of 
states analysis of Stampfl et ai^^. The character of the 
deeper lying bands is given in the caption of Fig. [TJ 

Focusing first on the band gaps presented in Table 
we note that both LDA and GGA underestimate all band 
gaps by more than 100%. The GGA band gaps are only 
marginally (~ 0.1 eV) larger than those of LDA. The 
reason is a combination of three factors: 1) in direct 
and inverse photoemission experiments electron addition 
and removal energies are probed, but the derivative dis- 



continuity of the exchange-correlation potential with re- 
spect to changes in the particle number is not taken into 
account in KS-DFT single-particle energy calculations— 
2) LDA and GGA are approximate exchange-correlation 
functionals, which 3) suffer from inherent self-interaction 
effects. The OEPx(cLDA) formalism also does not ful- 
fill criterion 2), but it is self-interaction free. This leads 
to a significant opening of the Kohn-Sham band gaps 
compared to those of LDA and GGA^i^SiiS as is ev- 
ident from Table HB Although the OEPx(cLDA) for- 
malism exhibits a derivative discontinuity^^'^^ and there- 
fore fulfills criterion 1) this is of no benefit in KS-DFT 
single-particle energy calculations. When the excitation 
energies are calculated by computing total energy differ- 
ences in OEPx(cLDA) between the N and the 7V±1 elec- 
tron system (frequently denoted A self-consistent field 
(ASCF) approach), the derivative discontinuity is taken 
into account properlySS,. In KS-DFT, however, the exci- 
tation energies are approximated by Kohn-Sham eigen- 
value differences of the A^-electron system alone. The 
derivative discontinuity does hence not enter the cal- 
culation and all states experience the same exchange- 
correlation potential. 

Having established that the removal of the self- 
interaction in the OEPx(cLDA)-KS approach is the dis- 
tinguishing feature compared to KS-LDA or KS-GGA 
calculations we will now illustrate how this leads to an 
opening of the band gap in ScN. For this it is illuminat- 
ing to start from the eigenvalues of the isolated Sc and 
N atoms, depicted in Fig. [31 The removal of the self- 
interaction in OEPx(cLDA) leads to a downward shift of 
all atomic states. Since the electrons in the 2nd shell of 
the nitrogen atom are more localized than the electrons 
populating the 3rd and 4th shell in scandium the self- 
interaction correction to the N 2p state is much larger 
than that of the Sc 3d state. Inspection of the differ- 



5 





-5 

^ -10 
> 

^-15 

at) 

i -20 
w 

-25 
-30 



: — — 4p 


... 4p : 


: ■■■ ~p 




L 27 

: Scandium Nitrogen 


3J : 
Gallium 



FIG. 3: Eigenvalue spectrum of the isolated Sc atom com- 
pared to N and Ga. Dotted lines show the levels in the LDA 
and solid lines in OEPx(cLDA). 



ence between the exchange potential in OEPx(cLDA) 
^^OEPx(cLDA)^ and LDA (w^dA) g^own in Fig. Ha) re- 
veals that the large relative shift of the atomic N 2p 
state also translates to the solid. Fig. (H^a) illustrates 
that ^o^P''('=LDA) significantly higher than v]^^^ in the 
Sc regions and lower around the N atoms. This difference 
in Vx leads to a significant charge density redistribution 
[shown in Fig. IHb)]. The charge transfer from the Sc 
to the N regions gives rise to an increase in the bond 
ionicity, which, in turn, leads to an opening of the band 
gap — consistent with our OEPx(cLDA) band structure 
calculations. 

In the II- VI compounds and group-III-nitrides this 
mechanism is also responsible for an opening of the band 
gap in OEPx(cLDA) compared to LDA, but it is com- 
plemented by a contribution arising from the coupling 
between the anion semicore d electrons and the 2p elec- 
trons of nitrogen. Taking GaN as example again the Ga 
?>d electrons of gallium are energetically lower than the 
2p electrons of nitrogen, while in Scandium the Sc 3d 
lie above the N 2p states (cf Fig. [3]). Not the anion 3c?, 
like in ScN, but the Ga 4s electrons thus form the lower 
conduction bands with the nitrogen 2p states. Since the 
OEPx(cLDA) shift of the N 2p states is larger than that 
of the Ga 4s states the bond ionicity and hence the band 
gap increase just like in ScN. In addition the Ga id elec- 
trons localize stronger in GaN when the self-interaction 
is removed in the OEPx(cLDA) approach^^. As a result 
the pd repulsion reduces and the valence bands are low- 
ered in energy leading to a further opening of the band 
gap^. 

We now turn to the quasiparticle band structure. In 
the OEPx(cLDA)-G'oW^o approach the band structure is 
calculated directly at the experimental lattice constant. 
However, the negative LDA band gap (see Table [Til im- 
pedes the application of this direct approach in the LDA- 
GqWq formalism with the gwst code, since in its current 
implementation^!^ a clear separation between conduc- 
tion and valence bands is required. Therefore, an indirect 



approach is adopted. First, LDA-GqWo calculations are 
performed at a lattice constant (oq = 4.75 A) larger than 
the experimental one, where the fundamental band gap 
in the LDA is small but positive. We then use the LDA 
volume deformation potentials (see Section [ill B[) to de- 
termine the corresponding LDA-GqW'o band gaps at the 
equilibrium lattice constant. Using the volume deforma- 
tion potentials of the LDA instead of the quasiparticle 
ones is a well justified approximation, as we will show in 
the next subsection. While this approach is in principle 
not limited to band gaps, it proves to be too cumbersome 
for a whole band structure calculation, because for every 
band structure point the corresponding deformation po- 
tential would have to be determined. 

The quasiparticle band structure is shown in Fig.[T]and 
the direct band gaps at the F and X-point are presented 
together with the indirect gap between F and X in Ta- 
ble [Til It is interesting to note that the LDA-GoM^o and 
OEPx(cLDA)-GoWo calculations, starting from the two 
extremes (negative band gap in LDA, 0.8 eV overestima- 
tion in OEPx(cLDA)), yield quasiparticle band gaps that 
agree to within 0.3 eV. Since the LDA-based calculations 
are close to the limit of metallic screening, whereas the 
OEPx(cLDA)-based calculations form the opposite ex- 
treme of starting from a completely self-interaction free 
exchange-correlation functional, we expect the results of 
a self-consistent GW calculation to fall in the range be- 
tween the LDA-GoWo and OEPx(cLDA)-GoWo calcula- 
tions. From these results we estimate the error bar asso- 
ciated with omitting self-consistency in GW to be of the 
order of 0.15 eV for ScN. Taking OEPx(cLDA)-GoWo 
results as lower and those of the LDA-GqWo as upper 
bounds, the arithmetic averages for Eg~^ , Eg~^ and 
E^~^ are 0.99, 3.62 and 2.02 eV, respectively. These 
are significantly lower than those from OEPx(cLDA) and 
the more approximate screened exchange calculations, as 
Tab . [Ill demonstrates . Our quasiparticle gaps clearly sup- 
port recent experimental findings of an indirect gap of 
0.9±0.1eV^^ and are at the lower bound of earlier mea- 
surements on samples with unintentionally high back- 
ground carrier concentration (1.3±0.3eV)l^. 

The fact that LDA-GqW^o and OEPx(cLDA)-GoWo 
calculations yield very similar quasiparticle band gaps is 
in disagreement with our previous observation for II- VI 
compounds and GaN^^. The difference between ScN and 
these compounds is that for the latter the cation semicore 
d shell is fully filled and the remaining s and p electrons in 
the semicore shell are much lower in energy. Taking GaN 
as an example the 3s electrons in the Gallium atom are 
approximately 100 eV and the 3p electrons approximately 
60 eV lower than in Scandium. Unlike in ScN the 3p de- 
rived bands therefore show no noticeable dispersion in 
GaN (cf Fig. [T]). Resolving these more strongly localized 
3s and 3p electrons in GaN with plane-waves will thus 
require significantly higher plane-wave cutoffs^^ than the 
80 Ry used in the present study for ScN. In a pseudopo- 
tential framework it would hence make sense to explic- 
itly include the d electrons of the cations in the II- VI 
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FIG. 4: The difference between (a) the OEPx(cLDA) and 
LDA exchange potentials (in Hartree) and (b) the electronic 
charge densities (in electrons/unit cell) of ScN for one of the 
square faces of the conventional rock-salt unit cell. Black 
circles denote Sc and white circles N atoms. 



compounds and group-III-nitrides as valence electrons, 
but to freeze the chemically inert semicore s and p elec- 
trons in the core of the pseudopotential. However, due to 
the large spatial overlap of the atomic semicore s and p 
with the d wavefunctions, core- valence exchange is large 
in these compounds. As a consequence core- valence ex- 
change is treated inconsistently when going from LDA to 
LDA-GqWo, if pseudopotentials are used in this fashion, 
because the exchange self-energy in the GW approach 



acts on the d electrons in the solid, but cannot act on 
the s and p electrons in the semicore shell, too. The re- 
sult is a severe underestimation of the LDA-GqWo band 
gaps and d-bands that are pushed energetically into the 
p-derived valence bands in the II- VI compounds^^i^i^. 

The only way to remedy this problem within LDA- 
GqWq is to free the electrons in question by perform- 
ing all-electron GqWq calculations'^^ or by using pseu- 
dopotentials that include the entire shell as valence 
electrons^i^^i^, which in the latter case introduces 
formidably high plane-wave cutoffs. If, on the other 
hand, OEPx(cLDA) is used for the ground state calcu- 
lation, then the exchange self-energy already acts on the 
semicore s and p states in the generation of the pseu- 
dopotential. Since the exchange self-energy can be lin- 
early decomposed into a core and a valence contribution 
no non-linear core corrections^ arise in the Hartree-Fock 
case and they are expected to be small for OEPx(cLDA) 
pseudopotentials'^^. We take the fact that the quasiparti- 
cle band structure in the OEPx(cLDA)-GoW^o approach 
agrees well with (inverse) photoemission data for these 
materials as indication that when switching from the lo- 
cal potential in OEPx(cLDA) to the non-local self-energy 
in OEPx(cLDA)-GoWo core- valence exchange is treated 
consistently, as long as OEPx(cLDA) pseudopotentials 
are used--. Since the semicore s and p states are less 
localized in ScN it is computationally feasible to include 
the entire 3rd shell of Sc as valence in the pseudopoten- 
tials (see Section |n| and thus to conduct a meaningful 
comparison between the LDA-GqWo and OEPx(cLDA)- 
GqWo calculations, which enables us to assess the error 
bar with respect to a self-consistent GW treatment. 

It remains to be added that in previous studies, where a 
significant starting point dependence in LDA-GqW^o com- 
pared to GGA-GqWo calculations was noted, this was ei- 
ther due to structural effects^ or significant differences 
in the ground state'^^ introduced when going from LDA 
to GGA. Since all calculations in this work were per- 
formed at the experimental equilibrium volume the KS 
band structures in LDA and GGA are very similar (cf. 
Tab. im and LDA-GqM^o and GGA-GqWo calculations 
yield essential the same result. 

Finally, Fig. [1] illustrates that the Sc 3p and Sc 3s 
bands are significantly lowered by the quasiparticle en- 
ergy calculations. We have argued recently^^ that this 
lowering can to a large degree be attributed to charge 
density relaxation effects arising from the removal of an 
electron from these states. These effects are accounted 
for in the GqWq approach, but not in DFT Kohn-Sham 
single-particle energy calculations and are larger for more 
localized states such as the Sc 3s and 3p bands. 



B. Effective masses and deformation potentials 

In this final part we present additional band structure 
parameters of ScN, namely the transverse (m^) and longi- 
tudinal (to* ) conduction band effective masses - at the X 
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1.92 
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0.126 


1.570 


1.36 


1.95 


2.03 


LDA^^ 


0.124 


1.441 









TABLE III: Transverse (m*) and longitudinal (m;*) effective 
masses of the conduction electrons at the X point (in units 
of mo) and different band gap volume deformation potentials 
{a'i) (in eV) for ScN. 




ln(V/V^) 

FIG. 5: Variation of the indirect band gap {E^~'^) and the 
direct band gap at the F-point (Sg of ScN with respect to 
the corresponding band gaps at the experimental equilibrium 
volume (Vb), as a function of ln(V/Vb). 

point - and the volume band gap deformation potentials, 
extracted from our quasiparticle energy calculations. For 
comparison we will also discuss the corresponding values 
obtained from LDA, GGA and OEPx(cLDA). 

The conduction band effective masses at the X point 
are calculated by fitting a quadratic function to the cor- 
responding band structure energies along the A (F-X) 
and Z (X-W) directions for m^* and m^, respectively. A 
small fc-point spacing of 0.01 in units of 27r/a yields con- 
verged effective masses, which are listed in Table Hill Our 
LDA results are in good agreement with those of Ref. WK 
As far as an experimental reference is concerned, we are 
only aware of one study, where a conduction band effec- 
tive mass between 0.1 and 0.2 mo has been reported^. 
Apart from the OEPx(cLDA) results, all our KS-DFT 
and quasiparticle energy calculations give a transverse 
effective mass in this range, while the longitudinal effec- 
tive mass is approximately one order of magnitude larger. 

To determine the volume deformation potentials 
for the band gaps Eg~^, Eg~^ and E^~^, the rele- 



vant band gaps are calculated at four different lattice 
constants around the experimental equilibrium one. For 
illustration, we show in Fig.[5]the variation of E^^^ and 
Eg~'^ , with respect to their values at the experimen- 
tal equilibrium volume, as a function of ln(V^/Vb). The 
corresponding band gaps are then fitted to a quadratic 
function of In V. The linear deformation potentials ob- 
tained in this way are listed in Table Hm To the best 
of our knowledge, no experimental or previous theoret- 
ical results are available for the deformation potentials 
of ScN. Based on our results we conclude that ay is al- 
most insensitive with respect to the computational ap- 
proaches considered in this article. The only exception 
is Uy^^ obtained using the OEPx(cLDA) approach. The 
fact, that the LDA curves in Fig. [5] are very close to the 
OEPx(cLDA)-G'oWo ones, a posteriori justifies the use 
of the LDA volume deformation potential for the indi- 
rect calculation of the LDA-GoIVo band gaps. 



IV. CONCLUSIONS 

Pseudopotential GqWq calculations based on Kohn- 
Sham density-functional theory calculations in both the 
LDA and OEPx(cLDA) have been performed for the elec- 
tronic structure of ScN in the thermodynamically stable 
rock-salt phase. To analyze the effects of exchange and 
correlation the atomic and electronic structure has been 
studied within DFT for several levels of approximations 
for the exchange-correlation functional (LDA, GGA and 
OEPx(cLDA)). In agreement with previous calculations 
for ScN, our LDA (OEPx(cLDA)) band gaps are under- 
estimated (overestimated) by about 100%. Despite this 
large difference, OEPx(cLDA)-GoWo and LDA-GqWo 
calculations for the quasiparticle band structure agree to 
within 0.3 eV. Our quasiparticle gap of 0.99±0.15eV sup- 
ports the recent observation that ScN has a much lower 
indirect band gap than previously thought. The main 
advantage of the OEPx(cLDA)-GoM^o approach lies in 
the fact that it facilitates a direct calculation of the elec- 
tronic structure of ScN at the experimental equilibrium 
volume, whereas for the LDA-GqWo calculation an indi- 
rect approach has to be taken due to the negative LDA 
band gap. 
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